A SEQUENCE OF MATRIX VALUED ORTHOGONAL 
POLYNOMIALS ASSOCIATED TO SPHERICAL FUNCTIONS 



I. PACHARONI AND P. ROMAN 

Abstract. The main purpose of this paper is to obtain an explicit 
expression of a family of matrix valued orthogonal polynomials {P„}„, 
with respect to a weight W, that are eigenfunctions of a second order 
differential operator D. The weight W and the differential operator D 
were found in [T2], using some aspects of the theory of the spherical 
functions associated to the complex projective spaces. We also find 
other second order differential operator E symmetric with respect to W 
and we describe the algebra generated by D and E. 



1. Introduction 

The theory of the harmonic analysis on homogeneous spaces is closely 
connected with the theory of special functions. This is apparent, for exam- 
ple, on the two dimensional sphere S'^ = SO(3)/SO(2), where the harmonic 
analysis with respect to the action of the orthogonal group is contained in 
the classical theory of the spherical harmonics. In spherical coordinates 
the spherical functions are the Legendre polynomials P„(cos0). Also the 
zonal spherical functions of the sphere S"" = SO(n + 1)/S0(n) are given, 

in spherical coordinates, in terms of Jacobi polynomials Pi"'"'*(cos^), with 
a = {n — 2)/2. More generally the zonal spherical functions on a Riemann- 
ian symmetric space of rank one can always be expressed in terms of the 
classical Gauss' hypergeometric functions, in the case of compact spaces we 
get Jacobi polynomials. 

As in the scalar case alluded above, in the matrix setting we also have 
these three ingredients: the theory of matrix valued spherical functions of 
any K-type, the matrix valued hypergeometric function and the theory of 
matrix valued orthogonal polynomials. In this paper we exhibit the interplay 
among these concepts in the case of the complex projective space -Pn(C) = 
SU(n + l)/U(n). 

The theory of matrix valued spherical functions goes back to [Mj and [5] , 
based on the foundational papers of Godement and Harish-Chandra. In [6] 
we find explicit expressions for spherical functions, of any K-type associated 
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to complex projective plane P2{C) = SU(3)/U(2). This is accomplished by 
associating to a spherical function <^ on G a vector valued function H defined 
on a complex affine plane C^, whose entries are given in terms of a special 
class of generalized hyper geometric functions p+i-Fp. 

The matrix valued hyper geometric function was studied in [15]. Let V be 
a d-dimensional complex vector space, and let A,B and C € End(y). The 
hypergeometric equation is 

(1) z{l - z)F"{z) + {C-z{A + B + I))F'{z) - ABF{z) = 0. 
If the eigenvalues of C are not in —No we define the function 

m=0 

where the symbol (C; A; B)m is defined inductively by 
(C;A;S)o = 1, 

(C; A- B)m+i = {C + m)-\A + m){B + m)(C; A; B)^, m>0. 

The function 2Fi[^^^;z^ is analytic on \z\ < 1 with values in End(l/). 
Moreover if Fq G F then F{z) = 2Fi[^^^ ; z^Fq is a solution of the hyper- 
geometric equation ([T]) such that F{0) = Fq. Conversely any solution F, 
analytic at z = is of this form. 

The theory of matrix valued orthogonal polynomials, without any con- 
sideration of differential equations goes back to pi)] and [11]. In [2], the 
study of the matrix valued orthogonal polynomials which are eigenfunctions 
of certain second order differential operators was started. The first explicit 
examples of such polynomials are given in [7j, [8] and [3]. 

Given a self adjoint positive definite matrix valued smooth weight function 
W = W{t) with finite moments, we can consider the skew symmetric bilinear 
form defined for any pair of square matrix valued polynomial functions P{t) 
and Q{t) by the numerical matrix 

(p,g) = / p{t)w{t)Q*{t)dt, 

jR 

where Q*{t) denotes the conjugate transpose of Q{t). This leads to the 
existence of a sequence of matrix valued orthogonal polynomials, that is a 
sequence where P„ is a polynomial of degree n with non singular 

leading coefficients and Pm) = if n ^ m. 

We also consider the skew symmetric bilinear form 

(2) {P,Q) = {P*,Q*r, 

and we say that a differential operator D is symmetric with respect to W 

if 



(3) {DP,Q) = {P,DQ), 

for all matrix valued polynomial functions P and Q. 
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Let D be an ordinary linear differential operator with matrix valued poly- 
nomial coefficients of degree less or equal to the order of derivation. If D 
is symmetric with respect to W then any orthogonal sequence {Pn}^ with 
respect to (■, •), satisfies 

(4) = P„*A„, 

for some numerical matrix A^. 

Assume that the weight function W = W{t) is supported in the interval 
(a, h) and let D he a, second order differential operator of the form 

(5) D = A2{t)^+Ai{t)j^+Ao{t), 

with matrix valued polynomial coefficients Aj{t) of degree less or equal to 
j. In [8] (see also [3j) it is proved that the condition of symmetry for D is 
equivalent to the following three differential equations 

AlW = WA2, 

(6) A\W = -WAi + 2{WA2)', 

A^W = WAq - {WAi)' + {WA2)", 
with the boundary conditions 

(7) YiraW{t)A2{t) = {) = Yiui{W{t)Ai{t)- A\{t)W{t)), ioi x = a,h. 

Finding explicit solutions of these equations is a highly non trivial task. 
In [3] and [1] the authors give some families of examples. In [12] one finds, 
for each dimension, a three parameter family of pairs {W, Z?} satisfying ([6]) 
and ([7]). These families arise from the representation theory of Lie groups. 
After the change of variable n = 1 — t, the main result in [12] reads: 

Theorem 1.1. Let a, (5 > -I, Q < k < 13 + I and i e n. Let D he the 

differential operator defined by 

with 

e e e 

2=0 i=l i=0 

e e-i 
V = J2i{a + p + i-k + l)Eu - -i){i + f3-k + l)^i,i+i. 

i=0 i=0 

Then the differential operator D is symmetric with respect to the weight 
matrix W{u) = (1 — u)"'u^ Z{u) given by 



ij=0 \r=0 
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Remark. Here, and in other parts of the paper, we use Eij to denote the 
matrix with entry equal 1 and otherwise. 

This theorem is obtained from the first few steps in the expUcit determina- 
tion of all matrix valued spherical functions associated to the n-dimensional 
projective space -Pn(C) = SU(n + 1)/U(n). The idea, also used in [6], is to 
cook up from a matrix valued spherical function a function H which depends 
on a single variable u. Using that the spherical functions are eigenfunctions 
of the Casimir operator of SU(n + 1) we deduce that, after an appropriate 
conjugation, H is an eigenfunction of an ordinary linear second order ma- 
trix valued differential operator D. The fact that this operator is symmetric 
with respect to the weight is a consequence of the fact that the Casimir 
operator is symmetric with respect to the L'^-inner product between matrix 
valued functions on SU(n -|- 1). At this point some readers may find useful 
to consult references [5], [Hj and [12]. 

One of the main purposes of this paper is to give explicit expressions 
of a sequence of orthogonal polynomials associated to the weight W given 
in Theorem II. li This is accomplished by studying the vector space V{X) 
of all vector valued polynomial solutions of the hypergeometric equation 
DF — XF = 0. This space is non trivial if and only if 

A = Xj{w) = -w{w + a + P + i + j + l)-j{a + p- k + l+j), 

for some if G No and j = 0,1, . . . , i. If the eigenvalues Xj{w) are all different 
then there exists a unique polynomial solution (up to scalars) of DF = XF. 
In Proposition 12.31 we compute, in the general case, the dimension of the 
space V{X). With this knowledge at hand, we construct a sequence of poly- 
nomials {Pw}, by choosing the j-th column of as a particular polynomial 
in V{Xj{w)). In Theorem 12.41 we prove that {Pw} is an orthogonal sequence 
of matrix valued polynomials such that DP^ = P^A^(Z)), where A„(L') is 
the real valued diagonal matrix 

0<j<i 

The matrix spherical functions associated to {G, K) = (SU(n -|- 1), U(n)) 
are eigenfunctions, not only of the Casimir operator, but also of any element 
in the algebra D{G)^ of all differential operators in G which are left and 
right invariant under multiplication by elements of G. In this case this 
algebra is a polynomial algebra in n algebraically independent generators, 
one of them can be taken to be the Casimir operator of G. For n = 2, in [6] 
the explicit expression of this set of generators was given and two differential 
operators D and E which commute were obtained. For a general n we do 
not have simple expressions for a complete set of generators of the algebra 
D{G)^ , beyond the Casimir operator. However in this paper we are able to 
find another second order differential operator E, which commutes with D 
and such that it is symmetric with respect to W (See Theorem 13. ip . The 
way in which we obtain this operator is different to the one used in [12] and 
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it is inspired in the operator E given in [13]. Here we only knew that such 
an operator should exist and after a trial and error process we find it and 
prove that it is symmetric. 

The sequence of matrix valued orthogonal polynomial constructed in The- 
orem [231 {-P«,} also satisfies EP* = P*A^{E) with 

e 

A^{E) = ^{-w{w + a + P + i + j + l){a - e + 3j) 

j=0 

-j{j + a + p-k + l)(a + 21 + Zk))Ejj. 

We also study the algebra generated by the differential operators D and 
E. In Theorem 15.31 we prove that it is isomorphic to the affine algebra of 
the following union of lines in C^: 

I 

J] (y - (q - £ + 3j> + 3j(^ -j + k){j + a + (5-k + l)). 

j=0 

Recently, in [T] this situation is considered in the case i = 2. The authors 
conjecture that the algebra generated by D and E coincide with the algebra 
of all differential operators that have the orthogonal polynomials Pw as 
simultaneous eigenfunctions. 

Acknowledgement. We would like to thank Prof. Juan Tirao for his 
continuous encouragement and for many useful comments and suggestions 
that helped us to improve this paper. 

2. Orthogonal polynomials associated to the pair {W, D} 

The aim of this section is to give explicitly a sequence of matrix valued 
orthogonal polynomials associated to the weight function W and the differ- 
ential operator D introduced in Theorem II. 1^ i.e. we construct a sequence 
{Pw} of orthogonal polynomials with respect to W, such that DP^ = P^^w, 
where A^(Z)) is a real diagonal matrix. 

The columns {Pii}j=o,...i of P^ are C^+^-valued polynomials such that 
DPi; = Xj{w)Pii and {Pi,,P^i) = 6w,w'^j,j'nw,j, for some positive real num- 
ber UwJ- 

2.1. Polynomial solutions of DF = XF. We start studying the C^^^- 
vector valued polynomial solutions of DF = XF. We will find all polynomials 
F{u) such that 

(8) n(l - u)F"{u) + (C - uU)F'{u) - {V + A)F(u) = 0, 

where the matrices C, U, V are given in Theorem 11.11 This equation is an 
instance of a hypergeometric differential equation studied in [H]. Since the 
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eigenvalues of C are not in — Nq the function F is characterized by Fq = F{0). 
For < 1 it is given by 

oo j 

(9) F{u) = 2Hi{^'l+^ ■,u)Fo = ^^[C;U;V + X],Fo, Fq G C^ 

1=0 

where the symbol [C;U;V + X]i is defined inductively by 
[C; U;V + X]o = l, 
[C; U;V + A]i+i = (C + i)-^ {{{U + i - 1) + V + X) [C;U;V + X]i, 
for all i > 0. 

There exists a polynomial solution of ([5]) if and only if the coefficient 
[C;U;V + X]i is singular for some i G Z. Let us assume that [C;U; V + X]w+i 
is singular and that [C;U;V + X]w is not singular. 

Since the matrix (C + w) is invertible, we have that [C;U;V + A]^+i is 
singular if and only if {w{U + w — 1) + V + X) is singular. The matrix 

= {w{U + w-l) + V + X) 

is upper triangular and 

(M^)jj = w{w + a + (3 + l + j + l)+i{a + (5-k + l+j) + X. 
Therefore [C]U]V + )^w+i is singular if and only if 

(10) X = Xj{w) = -w{w + a + (5 + (. + j + l)-j{a + (3-k + l + j), 
for some < j < 

We will distinguish the cases when the eigenvalues Xj{w) are all differ- 
ent (varying j or w) or when they are repeated. We start studying the 
polynomial solutions of ([8]) in the first case. 

Proposition 2.1. Assume that all eigenvalues Xj{w) are different. If X = 
Xj{w), for some j = 0,...,i, then there exists a unique Fq G C^^^ (up 
to scalars) such that F{u) = 2H1 (^'(7'*''^)''^) -^0 is a polynomial function. 
Moreover this polynomial is of degree w. 

Proof. We have already observed that for A = Xj (w) = —w{w + a + f3 + £ + 
j + 1) — j{a + (3 — k + 1 + j), the matrix [C, U,V + A]^+i is singular. Then 

the function F{u) = Yl'S=o irt^' ^ ^]iFo is a polynomial if and only if 
Fq is a vector such that 

(11) [C,U,V + XUFo£kei:{Mj; 

where = w{U + w - 1) + V + Xj{w). The matrix [C, U,V + A]^„ is 
invertible, hence Fq is univocally determined by an element in the kernel of 
Mu,. We have that 

(12) M^= {{i-j){a+l3-k+l+i+j+w)Eu-{e-i){P-k+l+i)Ei^i+i). 

o<i<e 
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Since all eigenvalues Xjiw) are different we have that 7^ Xj{w) — Xi{w) = 
{i — j){a + P — k + 1 + i + j + w) if i 7^ j, hence the dimension of the kernel 
of Mw is one. Explicitly (xq, xi, . . . , xe) G ker(Mu;) if and only if 

lY+iC^-i) ^ ii3-k+i+i)j_i forz-0 7 

|-]^3^ ^i-K ^) Ke-j) (a+/3+j+j+u,-fc+i),_>J loii-u, ...J, 

Xj+i = = ■ ■ ■ = Xi = 0, 

where we use {z)r = z{z + 1) . . . (z + r — 1), (z)o = 1. 

Hence, up to scalar, Fq is uniquely determined by (fTTj) and it is clear that 
F{u) = 2-ffi( ^ ' ; Fq is a polynomial of degree w with leading coefficient 
^[C, U,V + Xj{w)]wFo. This completes the proof of the proposition. □ 

Now we have to study the case when some eigenvalues are repeated, that 
is when there exist 10,10' G No and < j,j' < £ such that Xj{w) = Xji{w'). 
We start observing the following facts. 

Lemma 2.2. If Xj{w) = Xji{w') for some w,w' £ Nq and < < £ then 

i) We have w = w' if and only if j = j' . 

ii) If w' > w then j > j' + 1. 

Proof. If Xj{'w) = Xji{w') then 

{w' - w){a + [5 + £ + I + w + w' + ]) 

+(/ - + l3-k + l+j + j' + w) = ^. 

In particular if w' = w, we have (j' — j){a + 13 — k + 1 + j + j' + w) = 0. 
We observe that j / j' implies that a + fi — k + l+ j+j' + w > 0, because 
a > -1, /? - A; + 1 > 0, j + j' > 1 and tt; > 0. 

Similarly if j' = j we have {w' — w){a + (3 + i + l + w + w' +j) = 0. Since a > 
-1, P+£+l > and w+w'+j > 1 we obtain that {a+(3+£+l+w+w'+j) > 
and therefore w = w' . This completes the proof of i). 
For ii) we start from 

{w' -w){a + /3 + £ + l + w + w' = {j - j'){a + P -k + l + j + f + w), 

and we observe that the left hand side of this identity, as well as the factor 
(a + fi — k + l + j + j' + w) are positive numbers, by hypothesis, then we have 
j > j'. Finally suppose that j = j' + l then {w' — ■w){a + P + £ + w + w'+j) = 
{a + f3 — k + w + 2j), equivalently 

{w' - w - l)(a + (3 + £ + w + w' + j) = -{w' + £- j + k). 

The left hand side is non negative while the right hand side is negative 
because A; > 0, which is a contradiction. □ 

Let V{X) be the vector space of all C^^-'^-vector valued polynomials such 
that DP = XP. We observe that Proposition 12.11 said that if the eigen- 
values A = Xj{w) are all different the dimension of V{X) is one. The next 
proposition generalizes this result to the case when the eigenvalues Xj{w) 
are repeated. 
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Proposition 2.3. Let a,(5 > —1, < < /3 + 1 and let A = Xj{w), for 
some t(j G No- Then 

dim{P G V{X) : degP < w} 

= card {w' : < w' < w , X = Xj/^w'), for some < j' < £}. 

In particular 

dimV{X) = card {{w,j) : X = Xj{w)}. 

Proof. We have already observed that for A = Xj{w) the function F = F{u) 
is a polynomial solution of DF = XF if and only if F{u) = 2H1 (C, U, V+X)Fq 
with Fo G C^+^ such that [C, U,V + X]yjFo G ker(M^j), where 

Mw,j= Yl {ii-j)ia + P-k + l + i+j + w)E,i-{i-i){p-k + l + i)Ei,i+i) 
o<i<e 

We have that {i—j){a+P—k+l+i+j+w) 7^ if i 7^ j. Hence the dimension 
of ker(M^„j) is one. Moreover it is generated by (xq, • • • ,X£) G C^^-"^ such 
that 

= i-^r^ (H) for i = 0, . . . , - 1, 

(15) X, = 1 

Xj+i = Xj+2 = ■ ■ ■ = xe = 0, 

where we use {z)r = z{z + 1) . . . (z + r — 1), {z)q = 1. 

If the eigenvalue A is repeated s times and wi = m.u\{w G No : A = 
Xj{w),Q < j <i}, using Lemma [2^21 we can assume that 

A = Xj^{wi) = ■■■ = Xj^iws) 

with wi < W2 < ■ ■ ■ < Ws and ji > 32 + 1, J2 > is + 1, • • • , Js-i > is- 

For w = wi and j = ji the matrix [C, U,V + X]wi is invertible and 

Fq is univocally determined by an element in ker(M^u,^ j^), which is one 

dimensional, thus proving (jl4p in this case. 

Then to prove the proposition for any Wr we proceed by induction on 

1 < r < s. Thus let us assume that for 2 < r < s we know that 

{P G V{X) : degP < Wr^i} = r - 1. 

Let Mr = M^^j^. As we remarked 7^ P G V{X) is of degree Wr if and 
only if Po = P(0) 'satisfies ^ [C,U,V + A]^^Po G ker(Mr). 
Let 

[C7, [/, F + A]^, = NrMr-i . . . NiMiNo, 

where Ni are invertible matrices. The leading coefficient P^ of such a P is 
uniquely determined, up to scalar, by the condition 

MrNrMr^l . . . NiMiNoPq = 0, 

because we may assume that 

Pr = NrMr-l . . . NiMiNoPq = {xq, . . . , Xj^-i, 1,0,..., 0). 
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Now let us prove that there exists P G V^(A) of degree Wr, by constructing 
one by downward induction. 

Let Vr = {xq, . . . , a:j,,-i! 1) 0, . . . , 0) € ker(Mr) and let 6^ = N~^Vr. The 
equation br = Mr-iVr-i has a unique solution Vr-i of the form Vr-i = 
{zo, . . . , . . . , 0) because br = (yo, • • • , 0, . . . , 0) with yj^+i 7^ 

and Mr-i is upper triangular with a unique zero in the main diagonal in 
the jV-i-position. Similarly let br-i = N~_}-^^Vr-i, then there exists a unique 
Vr-2 = (^0) • • • ) tjr+'i, 0, ... 0) such that Mr-2Vr-2 = ^r-i- In this way we 
construct the sequence Vr, Vr-i, ■ ■ ■ ,vo such that 

Vr = Nrbr = NrMr-lVr-1 = NrMr-lNr-lMr-2Vr-2 = ■ ■ ■ 
= NrMr-l . . . NiMiNqVq 

Hence P = 2Hi(C, U,V + X)vr is a polynomial in V(\) of degree Wr- 
Now we observe that 

{P G V{X) : degP < Wr} = CP®{Pe V{\) : degP < Wr-i}. 

In fact it is clear that the right hand side is a direct sum contained in the 
left hand side. To prove the other inclusion we first observe that if i-* € ^(A) 
and dcgP < Wr then, as we saw, degP < Wr-i- If P € V{X) is of degree 
Wr then the leading coefficient of P is equal to the leading coefficient of tP 
for some t G C. Therefore P - tP e {P e y(A) : degP < Wr-i}. This 
completes the proof of the proposition. □ 

2.2. Matrix valued orthogonal polynomials associated to {W, D}. 

We want to construct a sequence {Pw}w>o of matrix valued orthogonal poly- 
nomials with respect to the weight function W, with degree of P^ equal to 
w, with non singular leading coefficient and that satisfies DP^ = P^A^, 
where Aw{D) is a real diagonal matrix. 

Then the columns {Pw}j=o,...e of P^ are C^'^^-valued polynomials such 
that Pi and P^, arc orthogonal to each other if {j,w) 7^ {j',w') and they 
satisfy that DP^ = Xj{w)Pw, where 

Xj{w) = -w{w + a + P + e + j + l)-j{j + a + l3-k + l), 

for each w G No, and j = 0, . . . , £. 

If an eigenvalue A = Xj{w) is not repeated, then we choose the unique 

Po G C^+i such that 



where denotes the i-th vector in the canonical basis of M + . Then we 
take 



(16) 




0<i<i 




00 



i=0 
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which is a polynomial function of degree w and satisfies 

DPi{u) = X^{w)Pi{u). 

(See Proposition 12. ip . 

If an eigenvalue A = Xj{w) is repeated we saw that, 

A = Xj,{wi) = Xj.,{w2) = ■■■ = Xj.iws), 

with wi < W2 < ■ • • < Ws and > > jV+i + 1, for 1 < r < s — 1. 

Let Vr = {P G ^(A) : degP < Wr}, for 1 < r < s. Then we saw, in the 
proof of Proposition 12.31 that 

/ Fi C ^2 C . • • C 

with dimVs = s Now we take, for each 1 < r < s 

/ P^;(n) = afl'i F^"- G Vr orthogonal to 

In this way, for each t(; G No we have defined i + 1 orthogonal polynomial 
functions P^, P^, . . . , P^ of degree w. 

Theorem 2.4. Let Pw{u) be the matrix whose rows are the vectors Pw{u). 
Then the sequence {-P«,(ti)}w,eNo an orthogonal sequence of matrix valued 
polynomials such that 

DPM = PUu)K, 
where = Yl^j=o Xj{w)Ejj . 

Proof. Let {w,j) / {w',j'). If Xj{w) ^ Xj:{w') then {Pi,P^,) = because 

D is symmetric. If Xj{w) = Xji{w') then {Pw,P^,) = by construction. 
Therefore the matrices P^, satisfies {Pw,Pw') = ii w ^ w' . 

On the other hand we have that for each w = 0, 1, 2 . . . the degree of 
Pw{u) is w and the leading coefficient of is the non singular triangular 
matrix 

-'"'"Z^V ^) \i^r) (a+/3+-r+s+u>-fe+l)r-s 
s<r 

This completes the proof of the theorem. □ 

3. The symmetry of the differential operator E 

The aim of this section is to exhibit another second order ordinary differ- 
ential operator which is symmetric with respect to the weight W. 

Theorem 3.1. Let a, (3 > -1, < k < (i + 1 and i e N. Let E be the 

differential operator defined by 

d^ d 
E={1- u){Qo + uQi)-^ + (Pq + uPi)-^ -ia + 2£ + 3k)V, 

with 

Qo = Yli=0 3^-E'i,i-l) 
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^0 = Eto ((« + + 1 + 2i) - 3k{i -i)- 3iiP -k + i))Eii 

- ELo ^(3^ + 3/3-3A; + 3 + ^ + 2a)^i,i_i, 
Pi = ELo -(a - ^ + 3i)(a + /? + ^ + z + 2)^^^ 

+ ELo 3(/3 - A; + 1 + i) - , 
y = Eto^(« + /? - A: + 1 + - Etd(^ - - A: + 1 + i)E,^i+i. 

Then E is symmetric with respect to the weight matrix W{u) = {l—u)'^u^Z[u), 
where Z[u) is given by 

^(-) = E f E (D (p n-n (i - uy-^u^^^) e,, 

ij=0 \r=0 J 

Proof. 

We need to prove that the equations ([6]) and ^ are satisfied. The equa- 
tions in ([6]) take the form 

(17) iQ*o + uQl)Z - Z{Qo + uQi) = 0, 

(Po* + uP^)Z + Z{Po + uPi) - 2Z{Qi -Qo- 2uZQi) 

^^^^ - 2(1 - n)Z'(Qo + uQi) - Mlz^b^2Z(go + uQi) = 0, 

+ (Po* + - Z'{Pq + uPi)- ZPi 

(19) + - T^) ((^0* + ^^Pi)Z - Z{Po + uP,)) 

- 2{a + 21 + 3k){ZV - V* Z) = 0. 

The ij-entry in the left hand side of the equation ()17p is 

u[a - £ + 3i)zij + 3{i + l)-Zi+ij - u{a - £ + 3j)zij + 3{j + l)zij+i = 0, 

because it is easy to verify that 

(i + l)zi+ij - {j + l)zij+i = u{j - i)zij. 

In order to prove the identity (jlSp we compute the ij-entry of the matrices 
involved there: 

{iP^+upm^J=^^^' E a)G)mrtn(i-^)'-^((^o),. 

r=max(i j) 

- {r - i){3i + 3p + e + 2a - 3k + 6) - {a - e + 3i){a + P + i + i + 2)'^ 

+ -'^^ E i'T) iT) CD Ct~r'7) (1 - uY-r 

r=max(i— 1 J — 1) 

(^(r - i + l)(3i + 3/3 + £ + 2a - 3A: + 6) + (a - ^ + 3i)(a + /? + ^ + i + 2))) 
+ u^^' E G-i) CT') ('^t.'"0 (1 - ^)'-''3(^ - . + l)(/3 + . - fc), 

r=max(i— l,j) 
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r 



(r - i)(3j + 3/3 + ^ + 2a - 3A: + 6) - (q - £ + 3j)(Q + /3 + £ + j + 2) 

E ri') (T) (SD (1 - 

r 

((r - i + l)(3i + 3/3 + £ + 2a - 3A; + 6) + (a - £ + 3i)(a + /3 + £ + j + 2)) 



{Z{Qi - Qo - 2uQi)) 



E (0 (p ri'') rtro (1 - + ^ - 3r) 



r 



+ E C't') d') (SO (T;'r) (1 - ^)'-^(3- + + 2a - 2£ + 3), 

r 

= E (I) (P ('^^") r t^P (1 - -Y-'-i^ - « - 3r) 

r 

+ E (IP (SiP ttrlr) (1 - -)^-^(« + /3)(3r + a - £ + 3), 



{{l-u)Z'{Qi+uQo)\j 

= n'^' E (P (P ('"TP ('i-;irP (1 - -)^-^^(- - - ^ + 3r) 



+ n'^' E (^IP (IP (SP (T;+r) (1 - -)^-^'- (3(- - J + 1)(^ - r + . + J-) 



+ (a-^ + 3j)(£-r + i + j-l)j. 
By using the previous results we get that the identity (fT8]l is equivalent 



to 

I 

E(r\ fr\ fl3+r-l\ (e+k-l^r\ _ ^ ,-,£-r 3{r+l){r+/3-k+l){e~r){2r+2^i-j) 
UAjA r A e-r A^ {r-i+l){r-j+l) 



+ (1 - (i) r (^+tj-p3(^ - j + k){j - ^) 



- E (^tP (IP (SP (T;'r) (1 - -)^-^3(^ - r + ^ - l)(2r + 2 - ^ - J) 

r=j-l 

= 0, 

which easily follows. 
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In order to prove the identity (fT9]) we compute 
{ZV - V*Z)ij 

= - ( - E (D (p r:r') c^tin (i - -)'-^(« +i-r+i) 

r 

+ E G) e^r') rt^o (i - + /? + ^ + j + ^ - r + 1)) , 

(p^z - zpo.- = i^ - jw^^-' ( E (0 (P ('""r') ('"-t^P 

r 

(1 - u)^"''+^(3i + 4a + 3/3 + 3A; + 6 - 3r + 5^ + 3j) 
- E (P (P rr') r '""P (1 - -)^-^(4a + M + 3k + 6- 3r)) , 

r 

. ^) an ^ _ ^^^^ ^ ^^^^^ ^ ^ , 

V — u) 

- z'{p, + np,))^ . = n^+^-i(. - ,) E (p (p r:r') rt.^-p (i - uy~^ 

(2(a + 2£ + 3A:)(q + ^ - r + 2) + 2q + ^ + 6 - 3A: - 3r) 

- .^+^--1 - J) E (D (p m n-n (i - -)^^^^^ ( 

r 

2(a + 2£ + 3/c)(a + /? + ^ - r + i + j + 3) + 3(i + j + /3 - 3A; + 2 - r - £)) . 

Now it is easy to verify that (fT9|) is satisfied. 

Finally the boundary conditions ([7]) can be easily check and this concludes 
the proof of the theorem. □ 

4. The algebra of differential operators 

Most of the results of this section are due to J. Tirao and they are taken 
from |16j. 

Let W = W{x) be a L X L matrix weight function with finite moments and 
let {Pn} be any sequence of matrix valued orthogonal polynomials associated 
to a weight function W. 

Let 

Vn = {Fe Ml^l{C)[x] : deg(F) < n} 

be the set of all matrix valued polynomials in the variable x of degree less 
or equal to n. 

Proposition 4.1. We have the following decomposition ofVn 

n 
j=0 
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Proof. It is clear that Yyj=o ^Lxl{C) is a subspace of Vn and that for 
n = they are the same. Let us denote by M„ the leading coefficient of P*. 

UH = ^„x"+A_ix"~^H h^o is a polynomial in K then H-P*M~^ An 

is a polynomial of degree < n — 1. Thus, by induction in n we obtain that 
H G Y^^=o PjMLxLiC). In order to prove that this sum is a direct sum 
we assume that PqAq + • • • + PnA^ = 0. By comparing, inductively the 
coefficients of x", x""^, ... x*^ we obtain that = • • • = = 0. □ 

Let "D be the algebra of all differential operators of the form 

(20) D = F,ix)£,+ Fs-i{x)£^ + • • • + ^i(^)^ + 
whith Fj a polynomial function of degree less or equal to j. 



Theorem 4.2. Let {Pn} be any sequence of matrix valued orthogonal poly- 
nomials associated to W . If D £ D is symmetric respect to W then DP* = 
Pn^n, for some matrix A„. 

Remark. We recall that D is symmetric with respect to W if {DP, Q) = 
{P,DQ), for all P,Q polynomials. The sequence {Pn} is orthogonal with 
respect to (,). The bilinear forms (,) and (,) are related by {P,Q) = 
{P*,Q*)*. 

Proof. Since D £ T> the operator D preserves the vector spaces Vn, for 
each n > 0. 

For n = we have that DPq G Vq, thus DPq = PqAq. By induction we 
assume that DP^ = Pj^j, for each < j < n — 1. By Proposition 14.11 we 
have that DP* = Y17=o ^i^i- Thus, for each < j < n — 1 we have 

{dp:,p*) = n=o(p:A^^p;) = Ei=i(^^,^)*^. = 

On the other hand, since D is symmetric we obtain 

{DP*n,P*) = {P:,DP*) = {P:,P*Aj) = ((P„,P,)A,)* = 0. 

Thus {Pj,Pj)*Aj = for each < j < n — 1, which implies that Aj = 
because the matrix {Pj,Pj) is non singular. Therefore DP* = Pn-^n and 
this concludes the proof. □ 

Given {Pn} any sequence of matrix valued orthogonal polynomials asso- 
ciated to the weight W, we define 
(21) 

V{W) = {D £V: DP* = P*An{D),yn > 0, for some matrix A„(Z))}. 

Proposition 4.3. We have 

(1) 'D{W) is a subalgebra of V which does not depend on the sequence 

{Pn}. 

(2) For each n G Nq, the function A„ : V{W) — > Mlxl{C) given by 
D I— > An{D) is a representation of the algebra V{W). 
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(3) The family {A„}„>o separates points of'D{W). That is, if Di and 
D2 are distinct points ofT>(W), then there exists no > such that 

K,{Di)^Kn,{D2). 

Proof. It is easy to verify that D{W) is a subalgebra of T>. To prove that 
it is independent of the sequence {Pn} we take other sequence of orthogonal 
polynomials {Qn}- Then Qn = AnPn, for some non singular matrix An- 
Then we have DQ* = DP*Al = P*An{D)A*^ = Q;T„(D), where T„(L>) = 

(^*)-ia„P)a;. 

If Di and D2 are in D{W) then 

D1D2P: = Di{P*K{D2)) = P„*A„pi)A„(Z)2). 

Hence A^pi^a) = K{Di)K{D2). 

Let us assume that there exists D € DiyV) such that A„(Z)) = for all 
n > 0. To prove (3) we have to verify that D = 0. By hypothesis we have 

that D = Yfi=QPi{x)£j satisfies DP* = 0, for ah n > 0. For n = we 
obtain FqP^ = 0, thus % = 0. 

By induction, we may assume that = for < i < j — 1, with j < s. 

Then = DP* = Y.i=i Piix)^-^^ = Fj{x)jlMj, where Mj is the leading 
coefficient of Pj, which is non singular. Therefore Fj =0. This concludes 
the proof. □ 

Corollary 4.4. The operators Di and D2 in the algebra T>{W) commute if 
and only if the matrices A„(Z)i) and A„(D2) commute for all n E Nq. 

Proof. By Proposition 14.31 (3) we have that D1D2 = D2D1 if and only 
if A„(i:>i_D2) = A„(i:>2-Di) for ah n. From Proposition (2) we get 
A„,{DiD2) = An{Di)AiD2). □ 

Proposition 4.5. Let {Qn} be the sequence of monic orthogonal polynomi- 
als. Let D = X]i=o-^«(^)^ such that DQ*^ = Q*n^n- Then 

(22) r„ = ^ [n]iJ^i for all n > 0, 

0<i<s 

where A\ is the coefficient of in the polynomial Fi. 

Remark. Here we are using the notation [n]j = n{n — \) . . . {n — i + 1) for 
n > 1, and [n]o = 1, for n > 0. 

Proof. From 

F,{u)^{u)=Q*n{u)Tn, 

0<i<s 

by comparing the monomials of degree n we get X]o<j<s["']«^i ~ ^n- D 

Remark 4.6. Observe that in particular, Proposition 14.51 implies that the 
eigenvalues is a polynomial function on n of degree less or equal to 
deg(D). 
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5. The operator E 

5.1. D and E commute. In this subsection we use the results described 
in Section [4] to give an elegant proof of the fact that the operators D and E 
commute. Of course we also can verify this by making the explicit compu- 
tations. 

Theorem 5.1. The differential operators D and E, introduced respectively 
in Theorems and \3.1\ commute. 

Proof. From Theorem 13.11 the operator E is symmetric with respect to 
the weight W . Thus E belongs to the algebra T>{W) defined in (pTj) (See 
Theorem 14. 2p . To see that D and E commute it is enough to verify that the 
corresponding eigenvalues commute. (See Corollary 14. 4p . 

Let {Qn} be the monic sequence of orthogonal polynomials. Then for 
any D S V{W), we have DQ"^ = Q'^VniD), where the eigenvalue Tn{D) is 
given explicitly in terms of the coefficients of the differential operator D (see 
Proposition 14. 5p . 

For the operators D and E introduced respectively in Theorems 11.11 and 
13.11 these eigenvalues are 

TniD) = -n{U + n-l)-V 

TniE) = -n{n - l)Qi + nPi - (a + 2£ + 3/c)y, 

where the matrices f7, V., Q\ and P\ are given in Theorems 11.11 and 13.11 
Explicitly we have 

r„(D) = -ELo("(^ + a + /3 + ^ + i + 1) + i(i + q + /? - A: + l))Eii 

+ Etd(^-^)(/? + ^-^ + i)^v+i 

^n{E) = -ELo("(« -i + ^i){n + a + p + £ + i + l) 

+(q + 2i + 3k)i{i + a + (3-k + l))Eii 
+ Ei=o(^ -i){P + i-k + l){a + 2£ + 3k + 3n)Si,i+i. 
Now it is easy to verify that 
(23) TniE) = {a + 2i + 3k + 3n)rniD) + 3n{e + k + n){n + a + l3 + i + l)L 

Thus the matrix Tn{E) commutes with Tn{D) and by Corollary 14.41 we have 
that D and E commute. □ 

5.2. The eigenfunctions of E. In Subsection 12.21 we give a sequence 
{Pw}w of matrix valued polynomials, which are orthogonal with respect 
to W and eigenfunctions of the differential operator D. The rows Pw of P^, 
are orthogonal polynomials of degree w and they satisfy DPi, = Xj{w)Pw. 

Since D and E commute, it follows that E preserves the eigenspaces of D. 
Therefore if an eigenvalue A = \j{w) has multiplicity one, then the vector 
valued polynomial Pw is also an eigenfunction of the differential operator E. 
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In the next theorem, we prove that this is true, even if the multiphcity of 
an eigenvalue is bigger than one. 

Theorem 5.2. The sequence {Pw}w of orthogonal polynomials associated 
to the pair {W, D} satisfies 

where Aw{E) = fij{w)Ejj, and 

^j{w) = -w{w + a+(3+£+j + l){a-£+3j)-j{j+a + (3-k + l){a + 2i+3k). 

Proof. Let {Q'!^}w>o be the sequence of monic orthogonal polynomials. 
Since E is symmetric with respect to the weight W, Theorem 14.21 savs that 
EQw = Qw^w{E) for some matrix ru,(i?). If = -P^^Jj, then we have 
that DP*Al, = P*Air^{D) and EP^Al, = P*Air^{E). Therefore 

A^D) = AlrUD){A*J-\ 

AUE)=AlrUE){A*J-'- 
Thus from ()23p we obtain that 

A^{E) = {a + 2i + 3k + 3w)A^{D) + 3w{£ + k + w){w + a + p + £ + 1)1. 

Observe that the fact that A^(L') is a diagonal matrix implies that A^(i?) 
is diagonal. Moreover the eigenvalue ^ij{w) = {Aw{E))jj is given by 

^j{w) = {a + 2l + 3k + 3w){-w{w + a + /3 + £ + j + l) 

-j{j + a + (3-k + 1)) + 3w{£ + k + w){w + a + 13 + £ + 1) 

= -w{w + a + [3 + £ + j + l){a - £ + 33) 

- j{j + a + 13 -k + l){a + 2£ + 3k). 

This concludes the proof of the theorem. □ 

5.3. The operator algebra generated by D and E. In this subsection 
we study the algebra generated by the differential operators D and E. 

Let C[x, y] be the algebra of all polynomials in the variables x and y with 
complex coefficients. 

Theorem 5.3. The algebra of differential operators generated by D and 
E is isomorphic to the quotient algebra C[x,y]/ (Q) , where (Q) denotes the 
ideal generated by the polynomial 

I 

Q{x,y) = J[{y-{a-£ + 3j)x + 3j{£ - j + k){j + a + f3-k + l)). 

j=0 

Proof. The algebra of differential operators generated by D and E is isomor- 
phic to the quotient algebra C[x, y]/I where I = {p € C[x, y] : p{D, E) = 0}. 

Since A^„ is a representation which separates points of ^{W) (Proposition 
|43|1 . we have that p{D, E) = if and only if 

AMD,E))=p{Ay,{D),AUE)) = 0, for all w. 
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Moreover, since the matrices A^(D) and A^(£^) are diagonal matrices, we 
have that p{Ayj{D),Aw{E)) = if and only if p{{Ayj{D))jj, {Aw{E))jj) = 
for all 0<j<L Thus the ideal I is 

I = {pe C[x, y] : piXj{w), fij{w)) = 0, for i = 0, 1, . . . i}. 

Let pj{x,y) be the polynomial 

Pj{x,y)=y-{a-i + 3j)x + 3j{i-j + k){j + a + P-k + l). 

It is easy to verify that pj{\j{w), iJ,j{w)) = 0, for all w > 0. Therefore 
Q{x,y) = Yl^j=QPj{x,y) belongs to the ideal /. 

On the other hand we have that any f € I vanishes in all points of the 
form (x, y) with y = (a — i + 3j)x + 3j{£ — j + k){j + a + (3 — k + 1), for 
each j = 0, . . . In fact if we let, 

a^=a-e + 3j bj = -3j{l-j + k){j + a + (3-k + l) {j = 0,1, . . . ,1) 

then we observe that the polynomial f{x, ajx + bj) has infinitely many roots, 
because f{\j{w),iij{w)) = and ^ijiw) = ajXj{w) + bj. 
Any polynomial in C[x,y] is also a polynomial in x and y — ax — b. Then 
it is clear that if p{x, y) = in the line y = ax + b then p is divisible by 
y — ax — b. 

Thus we have that if / belongs to the ideal / then / G rij^oipj) = iYljPj)- 
Therefore we have that the ideal / is generated by the polynomial Q{x,y), 
which concludes the proof of the Theorem. □ 
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